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We compare the groundstate energies of bosons and fermions with the same form of the Hamil- 
tonian. If both are noninteracting, the groundstate energy of bosons is always lower, owing to 
Bose-Einstein Condensation. However, the comparison is nontrivial when bosons do interact. We 
first prove that, when the hopping is unfrustrated (all the hopping amplitudes are non-negative), 
hard-core bosons still must have a lower groundstate energy than the fermions. When the hopping 
is frustrated, bosons can have higher groundstate energy than fermions. We prove rigorously that 
this inversion indeed occurs in several examples, using various techniques. 



Introduction. — Statistics of identical particles is one 
of the most fundamental concepts in quantum physics. 
One of the consequences of the particle statistics appears 
in the groundstate energy. For a system of free particles, 
the groundstate of bosons is obtained by putting all the 
particles in the lowest-energy state of the single-particle 
Hamiltonian, while the groundstate of fermions is ob- 
tained by filling the individual single-particle states up 
to the Fermi level (Pauli exclusion principle). Thus, the 
groundstate energy of bosons Eq and that of fermions 
for the same form of the Hamiltonian, satisfy 



(1) 



if the particles are non-interacting. 

The comparison becomes nontrivial when the particles 
do interact; the Bose-Einstein condensation is no longer 
perfect in interacting systems. Intuitively, it would be 
still natural to expect that Eq. (|T|) holds. However, re- 
cently an apparent counterexample was found numeri- 
cally 0, [ij ■ This motivates us to examine the fundamen- 
tal question how general the "natural" inequality ([T]) is 
and when it can be actually violated. 

To simplify the matter, in this paper we focus on the 
comparison of hard-core bosons with spinless fermions. 
The Hamiltonian is given as 



jkc]c k 



(2) 

where each site j belongs to a lattice A, rij = CjCj, and 
tjk = is assumed for j = k. The uniform (site inde- 
pendent) part of Hj is the chemical potential fj,. For a 
system of fermions, we identify Cj with the fcrmion an- 
nihilation operator fj satisfying the standard anticom- 
mutation relations. For a system of bosons, we identify 
Cj with the boson annihilation operator bj satisfying the 
standard commutation relations, with the hard-core con- 
straint rij = 0, 1 at each site. The hard-core constraint 
may be implemented by introducing the on-site interac- 
tion UJ2j n j( n j ~ 1) an d then taking U — > oo. 

We note that, the Hamiltonian (|2J) conserves the total 
particle number. Thus the groundstate can be defined 



for a given number of particles M (canonical ensemble), 
or for a given chemical potential /i (grand-canonical en- 
semble). The comparison between bosons and fermions 
can be made in either circumstance. 

Natural Inequality. - First we present a sufficient 
condition for the "natural" inequality (JTJ to hold. (See 
Ref. Q for a similar inequality for spinful fermions). 

Theorem 1. The inequality ([I]) holds for any given num- 
ber of particles M on a finite lattice A, if all the hopping 
amplitudes tjk are real and non-negative. 



Proof. Let us take the occupation number basis \<fi a ) = 
{{fij}}, where n° = M. The number operator rij 
has the same matrix element in this basis, for hard-core 
bosons and spinless fermions. On the other hand, for 
bosons, the matrix elements of each hopping term in the 
Hamiltonian ((2J is non-positive in this basis. Thus ma- 
trix elements for bosons H® b and those for fermions H^ b 
satisfy 



H 



at, 



-\K,\ 

njF 



(a^b) 
(a = b) 



(3) 



Let |*o)f = S a ^ a l^ a ) F b e the (normalized) ground- 
state for fermions. To obtain an upper bound on the 
groundstate energy of bosons, we consider the trial state 
|T'o)b = J2 a l^aW^B, where \4> a )B is the basis state for 
bosons corresponding to |</> a )F- Then Eq. ((3]) implies 

Ef < b(*o|H|* )b < f(*o|H|* )f = El. (4) 

□ 

As a simple corollary, the groundstate energies for a 
given chemical potential fi also satisfy Eq. ([T]). 

The sign of hopping amplitudes in a system of bosons 
is related to frustration. This is clear if one maps the sys- 
tem of hard-core bosons to the quantum spin system with 
S = 1/2; positive hopping amplitudes tjk correspond to 
ferromagnetic XY interaction. Theorem Q] means that, 
if there is no frustration in the hopping term, hard- 
core bosons always have lower energy than corresponding 
fermions. This also gives an alternative understanding 



of why the fermions have higher energy in the "unfrus- 
trated" case: statistical phases of fermions introduce a 
kind of frustration to the system. 

Let us now discuss how the natural inequality ^ can 
be violated. According to Theorem [TJ in order to realize 
the violation, it is necessary to introduce a frustration 
by choosing some amplitudes tjk to be negative or com- 
plex. While the presence of frustration is not a sufficient 
condition, we will demonstrate that the violation indeed 
occurs in several concrete examples. In the following, for 
simplicity we only discuss tight-binding models of hard- 
core bosons and corresponding nonintcracting fermions, 
setting Vjk = 0. 

Particles on a Ring. — We begin with a simple but in- 
structive example in one dimension: tight-binding model 
on a ring H = — X)jli( c j c j+i + H.c). The hard-core 
boson version of this model, which is equivalent to the 
S =1/2 XY chain, can be mapped to the model of non- 
interacting fermions on the ring by the Jordan- Wigncr 
transformation 0, 0|- Thus the hard-core bosons and 
fermions are almost equivalent in this case. 

However, a care should be taken on the boundary con- 
dition when we discuss the ring of finite length. For 
simplicity, here we assume that the number of sites 
N is an integral multiple of 4, and the number of 
particles M = N/2 (an even integer by assumption). 
Then the hard-core bosons with the periodic (antipcri- 
odic) boundary condition is mapped to nonintcracting 
fermions with the antiperiodic (periodic, respectively) 
boundary conditions. Nonintcracting fermions on a ring 
have a groundstate energy density lower by 0(1/N 2 ) for 
the antiperiodic boundary condition, compared to the 
periodic boundary condition jjsj]. This implies that the 
hard-core bosons have lower energy than fermions on 
a ring with the periodic boundary condition, conform- 
ing to Theorem Q] since all the hopping amplitudes are 
non-negative. On the other hand, the same result im- 
plies that, under the antiperiodic boundary condition, 
the hard-core bosons have higher energy than fermions. 
Imposing the antiperiodic boundary condition is equiv- 
alent to introducing a 7r-flux inside the ring, which can 
cancel the effect of the statistical phase so that the in- 
equality ([1]) is indeed inverted. The energy difference on 
the ring, however, vanishes asymptotically in the thermo- 
dynamic limit N — > oo. Thus, wc shall seek for different 
examples where the hard-core bosons have higher energy 
than fermions in the thermodynamic limit. 

Two Dimensional System with Flux. — A natural sys- 
tem to consider would be a two-dimensional lattice with 
flux. Under the periodic boundary condition, the total 
flux is quantized as integer flux quanta (the unit flux 
quantum is 2tt in our unit). For the finite square lattice 
with N plaquettes, the flux per plaquette $ is quantized 
in unit of 2ir/N . Such a uniform flux can be represented 
using the string gauge 0- We obtained the groundstate 
energy of hard-core bosons and fermions with various 




Figure 1. Difference of groundstate energy density between 
the hard-core bosons and fermions (Ae = E$/N - Eq /N) 
on the \/26 x \/26 square lattice with $ flux per plaquette 
and n e particle per site. The natural inequality (ffj holds in 
white region, while its violation is color coded. Statistical 
transmutation is expected along the two solid lines $ = n e 
and $ = 1 — n e . 

densities of particles and various values of flux using exact 
numerical diagonalization, for square lattices with peri- 
odic boundary conditions, up to 26 sites. The result for 
the y/26 x -s/26 square lattice is shown in Fig. [TJ 

We find that the "natural" inequality ((T|) is violated 
in a region of the phase diagram. In particular, the in- 
version is significant along the diagonal lines $ = n e 
and $ = 1 — n e - These lines are precisely where the 
statistical transmutation between the hard-core boson 
and the fermion is expected to occur 0, Q . Namely, in 
the mean-field level, one flux quantum can be attached 
to each particle, transforming fermions into bosons and 
vice versa, at the same time eliminating the background 
field. At zero field, the frustration is absent and hard- 
core bosons have lower energy than fermions. Thus, the 
statistical transmutation implies that, hard-core bosons 
have higher energy than fermions on the lines $ = n e 
and $ = 1 — n e . While this argument is not rigorous, our 
numerical result supports the statistical transmutation 
scenario. (For a related discussion for spinful electrons, 
see Ref. [l0(.) Numerical results for the square lattices of 
various sizes up to 26 sites (not shown) suggest that the 
energy difference is nonvanishing in the thermodynamic 
limit. Furthermore, in the following, we rigorously prove 
the inversion of the groundstate energies for several lat- 
tices in the thermodynamic limit. 

Flat Band. — Some models of noninteracting fermions 
have completely flat bands as the lowest energy band. In 
these models, as long as all the fermions occupy the low- 
est flat band, Pauli exclusion principle plays no role in 
increasing the groundstate energy. Thus, such flat band 
models would have better chance to realize the inversion 
of the groundstate energies. Here we show that the in- 
equality (TTJ) is indeed violated in a few examples with flat 
bands. 

First we discuss the delta-chain model, for which the 
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violation of Eq. ([I]) was numerically found for small clus- 
ters P, Q • The Hamiltonian of the model can be written 
in the following form 11 , TiJ : 



H 



N 

0=1 



(5) 



where the a-operator, which acts on each triangle, is de- 
fined as Oj = C2j_i + v / 2c 2 _, + C2j+i. Here we use a 
periodic boundary condition to identify C2N+1 with c\. 
The Hamiltonian % corresponds to a model with nega- 
tive hopping amplitudes tjk (as defined in Eq. ((5J)), which 
lead to frustration. 

The model in the single-particle sector has two bands. 
The lower flat band with zero energy is spanned by states 
annihilated by a/s. We note that the Hamiltonian ([5]) 
is modified from that in Ref. [U by a constant chemical 
potential, so that the flat band has exactly zero energy. 
Thus the groundstate energy of the fcrmionic version of 
the model ([3]) is zero as long as the filling fraction (par- 
ticle number per site) v satisfies v < 1/2. On the other 
hand, the groundstate energy Eq 3 of bosons is zero as 
long as v < 1/4 since the localized zero-energ y st ates do 
not overlap each other in this range of filling [ 1 31 ] . 

Now let us derive a nontrivial lower bound for Eq for 
filling fractions v > 1/4. We decompose the model into 
clusters, each containing p unit cells: 



H 



N/p-l 

£ 

71 = 



N/p 

y ] a np a n 



(6) 



where % 



(/.<) 



T p ~ x a} 



np +jO>np+j is the Hamiltonian 



for the (n + l)-th cluster. Since the second term 

J2n= P i a np a npi which connects neighboring clusters, is 

positive semidefinite, the groundstate energy Eq 3 of the 

first term H = YZ=i l ^ satisfies Eq < Eq. U is a 

(p) 

sum of mutually commuting cluster Hamiltonians TL„ . 
Thus Eq 3 is simply given by the sum of groundstate en- 
ergies of each cluster. The particle number within each 
cluster is also conserved separately in H. Let us choose 
p = 4 as in Fig. [5J The cluster contains 8 sites and the 
number of particles is to = 0, 1, . . . , 8. We obtained the 
groundstate energy Cq (to) of the cluster in each par- 
ticle number to sector, by numerical exact diagonaliza- 
tion. It is found that eg (m) = for < m < 3, and 

4 4) (m) > A^c = 0.372605 for 4 < to < 8. 

If we consider the filling fraction in the range 3/8 < 
v < 1/2, it follows from Dirichlet's box principle that 
there is at least one cluster which contains 4 or more 
particles. Thus, in this range, Eq > A^C for any system 
size N, while Eq 7 = 0. Therefore, the inversion of the 
groundstate energies holds also in the thermodynamic 
limit. 

This method can be easily extended to other lat- 
tices. For example, the standard nearest-neighbor hop- 





Figure 2. Decomposition of the delta-chain Hamiltonian to 
clusters. Each cluster contains p — 4 unit cells (8 sites), in- 
cluding one decoupled site at the top of the dashed triangle. 
The hoppings on dashed triangles 5Zn=i a np a np are omitted 
from the Hamiltonian. Since the omitted terms are positive 
semidefinite, the groundstate energy for the decoupled Hamil- 
tonian H gives a lower bound for the groundstate energy of 
the original Hamiltonian T-i. 



ping model on a kagome lattice is described by 

n = E a L a Ao +E a V Q a V Q ' ( ? ) 

a a 

where ^ a and Y a are elementary triangles pointing up 
and down respectively, of the kagome lattice. We de- 
fine = c Ql + c Q2 + c a3 , where 01,2,3 refer to the 
three sites belonging to fa a , and likewise for ay . The 
fermionic version of the model has three bands, the low- 
est of which is a flat band at zero energy EMU. Thus 



Eq = when v < 1/3. For the groundstate energy of 
the bosonic version, we can use the cluster decomposi- 
tion technique similar to what we have discussed above 
for the delta-chain. Let us choose the 12-site cluster of 
the "Star of David" shape. The groundstate energy eo of 
each cluster is zero in the fixed particle number sectors 
with to < 3, but is positive in the sectors with to > 4. 
Thus, invoking Dirichlet's box principle again, Eq. (TTJ) is 
violated for filling fraction 1/4 < v < 1/3. 

The outcome of the above argument depends on the 
cluster size taken. In fact, the range of filling fraction 
v for which we have proved the violation of Eq. (fTJ) is 
not optimal. In an alternative approach generalizing the 
techniques used in the context of flat-band ferromag- 
netism [ll[ and of frustrated antiferromagnets near the 
saturation field [13j , w e can extend the region of v where 
Eq. (fTJ) is violated [Iff. For example, in the delta-chain 
model, Eq. ([Tj is violated for 1/4 < v < 1/2. 

Anderson's argument. — In order to deal with more 
general models, we develop an alternative approach. Let 
us discuss the square lattice with 7r-flux per plaquette, at 
zero chemical potential. The Hamiltonian reads 



<0,k> 



(8) 



where tjk = ±1 as specified in Fig. [31 We note that Lieb 
has shown that 7r-flux minimizes the groundstate energy 
of fermions on the square lattice [ItJ ■ On the other hand, 
an argument similar to the Proof of Theorem [1] can be 
used to prove a lattice version of Simon's theorem on 
diamagnetism of bosons [3]. Namely, for bosons, intro- 
duction of flux always increases the groundstate energy. 
These, together with the statistical transmutation argu- 
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Figure 3. The square lattice with n flux in each plaquette. 
We choose the gauge so that the hopping amplitude tjk is 
+ 1 on black links, and —1 on blue ones. The brown cross 
represents a cluster of 12 sites. 



merit discussed earlier, suggest a possibility of violation 
of Eq. (fT|) with 7r-flux per plaquette. 

The fermionic version of this model does not have a 
flat band. Nevertheless, from the exact dispersion re- 
lation, the groundstate energy at half filling (fi = 0) is 
obtained exactly as Eq ~ — 0.958091A in the limit of 
large A. Now we turn to the "grandcanonical" ground- 
state energy of the corresponding boson model at the 
same chemical potential (fi = 0). Here we use Ander- 
son's argument 191-21 by writing the Hamiltonian as a 
sum of sub-Hamiltonians: 



U = ^Th a . 



(9) 



The groundstate energy Eq of H. satisfies Eq > e o> 
where £q is the groundstate energy of h a . For the present 
problem, we can choose h a as 



H.c.), 



(10) 



where + Q refers to a cross-shaped cluster of 12 sites as 
shown in Fig. [3] Wc consider all the clusters with the 
same pattern of hopping amplitudes within the cluster, 
in the square lattice. As a consequence, each cluster over- 
laps with 4 other clusters and each link appears in two 
different clusters when periodic boundary conditions are 
imposed. The factor 1/2 in Eq. (|10[) compensates this 
double counting. The grandcanonical groundstate energy 
of the cross-shaped cluster is obtained by exact diagonal- 
ization as Cq = —3.609035. 

Since there are A/4 such clusters in the square lattice 
of N sites, we obtain 

#o /N > -3.609035/4 = -0.902259 > Eq/N. (11) 

Thus the inversion of the groundstate energies is proved 
for the square-lattice 7r-flux model with /i = 0. This 
relation, expected from the statistical transmutation ar- 
gument discussed earlier, is now proved rigorously. 



This argument is not restricted to two dimensional 
systems. Let us consider the standard tight-binding 
model on a pyrochlore lattice: % — k)( c j c k + H.c), 
where (j, k) runs over all pairs of nearest-neighbor sites 
in the three dimensional pyrochlore lattice. Again we 
set the chemical potential fi — 0. We note that this 
model has frustrated hoppings with this choice of the 
sign. The model in the single-particle sector has two 
degenerate flat bands at the energy £ = — 2 and two 
dispersive bands touching the flat bands [Hj]. Thus 
for fermions, the groundstate energy at fi = satisfies 
El < -2(N/2) = -A, where A is the number of sites of 
the lattice. 

The hard-core boson version of this model can be de- 
composed as Eq. (j9]) with h a = (1/4) J2{j,k)eTD a ( c j c k + 
H.c), where TD Q refers to each dimer of elementary 
tetrahedra of the pyrochlore lattice sharing a vertex 
(site) . Here we count dimers in any direction; each tetra- 
hedron (and thus each link) belongs to 4 dimers. The 
factor 1/4 in the definition of h a is introduced to com- 
pensate the overcounting. The groundstate energy of a 
tetrahedra dimer is obtained by exact diagonalization as 
e C Q l = -(2+\/2)/4 = -0.853554. Since there are A dimers 
of tetrahedra, the groundstate energy of bosons at // = 
satisfies 



Eq/N > -(2 + \/2)/4 > El/N. 



(12) 



Thus we have proved the violation of Eq. ([T]) for the 
simple tight-binding model on the three-dimensional py- 
rochlore lattice. 

Conclusions. — We have investigated the fundamen- 
tal question whether the groundstate energy of hard-core 
bosons is lower than that of fermions on the same lattice. 
We have proved that the former is indeed lower than the 
latter, as naturally expected, if there is no frustration in 
the hopping. However, the inequality can be reversed in 
the presence of frustration, and we have demonstrated 
that it is actually the case in several concrete examples. 
There is a close connection of the present problem to 
many apparently unrelated concepts in quantum many- 
body problem, including Simon's universal diamagnetism 
of bosons, Lieb's optimal flux for fermions, statistical 
transmutation, flat-band fcrromagnetism, and antiferro- 
magnets near the saturation field. More details of our 
analysis will be presented in a separate publication fl6| . 
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